Lorentz and CPT violation in the neutrino sector 
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We consider neutrino oscillations in the minimal Standard-Model Extension describing general 
Lorentz and CPT violation. Among the models without neutrino mass differences is one with two 
degrees of freedom that reproduces most major observed features of neutrino behavior. 
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Quantum physics and gravity are believed to combine 
at the Planck scale, mp ~ 10^^ GeV. Experimentation 
at this high energy is impractical, but existing technol- 
ogy could detect suppressed effects from the Planck scale, 
such as violations of relativity through Lorentz or CPT 
breaking [1]. At experimentally accessible energies, sig- 
nals for Lorentz and CPT violation are described by the 
Standard-Model Extension (SME) [2], an effective quan- 
tum field theory based on the Standard Model of parti- 
cle physics. The SME incorporates general coordinate- 
independent Lorentz violation. 

The character of the many experiments designed to 
study neutrino oscillations [3] makes them well suited 
for tests of Lorentz and CPT symmetry. The effects of 
Lorentz violation on propagation in the vacuum can be- 
come more pronounced for light particles, and so small 
effects may become observable for large baselines. Ap- 
plying this idea to photons has led to the best current 
sensitivity on any type of relativity violation [4] . 

In this work, we study the general neutrino the- 
ory given by the minimal renormalizable SME [2]. In 
this setup, as in the usual minimal Standard Model, 
SU(3)xSU(2)xU(l) symmetry is preserved, the right- 
handed neutrino fields decouple and so are unobservable, 
and there are no neutrino mass differences. The neutrino 
behavior is contained in the terms 

C D jiLal^ La - {aL)^abLa7''Lb 

+ ^t{cL),.uabLal'' D"" Lb , (1) 

where the first term is the usual Standard-Model kinetic 
term for the left-handed doublets La, with index a rang- 
ing over the three generations e, ^, t. The coefficients 
for Lorentz violation are {aL)ij,ab, which has mass dimen- 
sion one and controls the CPT violation, and {cL)^iiab^ 
which is dimensionless. It is attractive to view these co- 
efficients as arising from spontaneous violation in a more 
fundamental theory [5] , but other origins are possible [1] . 

The Lorentz- violating terms in Eq. (1) modify both in- 
teractions and propagation of neutrinos. Any interaction 
effects are expected to be tiny and well beyond exist- 
ing sensitivities. In contrast, propagation effects can be 
substantial if the neutrinos travel large distances. The 
time evolution of neutrino states is controlled as usual 
by the effective hamiltonian (/ioff)ab extracted from Eq. 
(1). The construction of {hes)ab is complicated by the un- 
conventional time-derivative term but can be performed 



following the procedure in Ref. [6]. We find 

{heff)ab = \p\Sab + ^ [(aL)>^ " (clT" Pf,Pt.]ab- (2) 

To leading order, the 4- momentum is = (|p|; —pj. 

The analysis of neutrino mixing proceeds along the 
usual lines. We diagonalize (ft.cff)afc with a 3 x 3 unitary 
matrix Ucs, /icff = C^cff^-cff C^ff, where E^s is a 3 x 3 diag- 
onal matrix. There are therefore two energy-dependent 
eigenvalue differences and hence two independent oscil- 
lation lengths, as usual. The time evolution operator is 

Svc^b^^) — {UlffS~'^^°"*Uctf)ab, and the probabihty for a 
neutrino of type b to oscillate into a neutrino of type a 
in time t is P^^^^^{t) = \S^^^^{t)\'^ . 

The CPT-conjugate hamiltonian h^g"^ is obtained by 
changing the sign of ql- Under CPT, the transition am- 
plitudes transform as Sy^^^^{t) ^ Sl^fj^{-~t), so CPT in- 
variance implies Pu^^^^{t) — Pj^^^p^{t). Note that the 
converse is false in general [7]. For instance, the model 
described below violates CPT but satisfies the equality. 

Since oscillations are insensitive to terms proportional 
to the identity, each coefficient for Lorentz violation in- 
troduces two independent eigenvalue differences, three 
mixing angles, and three phases. The minimal SME 
(without neutrino masses) therefore contains a maximum 
of 160 gauge-invariant degrees of freedom describing neu- 
trino oscillations [8]. Of these, 16 are rotationally invari- 
ant. The existing literature concerns almost exclusively 
the rotationally invariant case [9-12], usually with either 
o-L or cl neglected and in a two-generation model with 
nonzero neutrino masses. A wealth of effects in the gen- 
eral case remains to be explored. 

The presence of Lorentz violation introduces some 
novel features not present in the usual massive-neutrino 
case. One is an unusual energy dependence, which can be 
traced to the dimensionality of the coefficients for Lorentz 
violation. In the conventional case with mass-squared 
differences Ato^, neutrino oscillations are controlled by 
the dimensionless combination Am'^L/ E involving base- 
line distance L and energy E. In contrast, Eq. (2) shows 
that oscillations due to coefficients of type and are 
controlled by the dimensionless combinations a^L and 
clLE, respectively. 

Another unconventional feature is direction-dependent 
dynamics, which is a consequence of rotational- symmetry 
violation. For terrestrial experiments, the direction de- 
pendence introduces sidereal variations in various ob- 
servables at multiples of the Earth's sidereal frequency 
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~ 27r/(23 h 56 m). For solar- neutrino experiments, it 
may yield annual variations because the propagation di- 
rection differs as the Earth orbits the Sun. Both types of 
variations offer a unique signal of Lorentz violation with 
interesting attainable sensitivities. For solar neutrinos 
LE ~ 10^^, so a detailed analysis of existing data along 
the lines of Refs. [14] might achieve sensitivities as low as 
10~^^ GeV on ol and 10~^^ on cl in certain models with 
Lorentz violation. These sensitivities would be compara- 
ble to the best existing ones in other sectors of the SME 
[4, 15-21]. 

The coefficients for Lorentz violation can also lead to 
novel resonances, in analogy to the MSW resonance [22]. 
Unlike the usual case, however, these Lorentz-violating 
resonances can occur also in the vacuum and may have di- 
rectional dependence [23] . Note that conventional matter 
effects can readily be handled within our formalism (2) by 
adding the effective contributions (aL,eff)ee = Gpi'^rie — 
nn)/V^ and (aL,cff)°^ = {aL.cs)rT = -Gprin/V^, where 
ne and n„ are the number densities of electrons and neu- 
trons. The contributions to hes from matter range from 
about 10-2° GeV to IQ-^^ GeV. This range is within 
the region expected for Planck-scale Lorentz violation, 
so matter effects can play a crucial role in the analysis. 

An interesting question is whether the introduction of 
Lorentz violation may help explain the small LSND ex- 
cess of De [24] . Usually, two mass-squared differences are 
invoked to explain the observations in solar and atmo- 
spheric neutrinos, but LSND lies well outside the region 
of limiting sensitivity to these effects. Possible solutions 
to this puzzle may arise from the unusual energy and di- 
rectional dependences of Lorentz violation. An explana- 
tion of LSND requires a mass-squared difference of about 
10-19 GeV2 = IQ-i eV2, an of about 10'^^ GeV, or 
a Cl of about IQ-^^. Any of these would affect other 
experiments to some degree, including the MiniBooNE 
experiment [25] designed to test the LSND result. 

To illustrate some of the possible behavior allowed by 
the SME, we consider a two-coefficient three-generation 
case without any mass-squared differences, but incorpo- 
rating an isotropic cl with nonzero element ^{cL)e^ = 
2c and an anisotropic gl with degenerate nonzero real el- 
ements (aL)f^ = {aL)eT = ajypl. The coefficients are un- 
derstood to be specified in the conventional Sun-centered 
celestial equatorial frame (T, X, y,Z), which has Z axis 
along the Earth rotation axis and X axis towards the 
vernal equinox [13]. In what follows, we show that this 
simple modcil. which we call the 'bic;yclc' model, suffices 
to reproduce the major features of the known neutrino 
behavior other than the LSND anomaly, despite having 
only two degrees of freedom rather than the four degrees 
of freedom used in the standard description with mass. 

Diagonalizing the hamiltonian for the model yields 

P,^^,^ = l-4sin2 0cos2 0sin2(A3iL/2), 
P^^^^^ = Py,^u^ =2sin2 6'cos2 6»sin2(A3iL/2), 



P,^^,^ = P,^^,^ = l-sin2 0sin2(A2iL/2) 

- sin^ e CQf? e sin2(A3iL/2) 

- cos2 6'sin2(A32L/2), 
P^^«^^ = sin2 6isin2(A2iL/2) 

- sin^ e cos^ e sin^ ( A31 L/2) 

-Fcos2 6»sin2(A32L/2), (3) 

where 

A21 = VicE)^ + (dcosQ)^ +cE, 
A31 = 2v'(c^)2 + (acose)2, 
A32 = ^/ {cE)^ + {a cos ey-cE, 
sin'^e = i[l - cEI^{cEf (acose)2], (4) 

and where we define the propagation direction by the 
unit vector j5 = (sin cos $, sin sin $, cos 0) in polar 
coordinates in the standard Sun-centered frame. These 
probabilities also hold for antineutrinos. 

The qualitative features of the model can be under- 
stood as follows. At low energies, a causes oscillation of 
into an equal mixture of and v^. At high ener- 
gies, c dominates and prevents mixing. For definite- 
ness, we take c > 0. At energies well above the criti- 
cal energy = |a|/c, svc? d vanishes and the probabili- 
ties reduce to a maximal-mixing two-generation v^, ^ Vt 
case with transition probability P^^^y^ ~ sin2(A32i/2), 
A32 — 0? cos^ Q/2cE. The energy dependence in this 
limit is therefore that of a conventional mass-squared 
difference of Arn|, = o^cos^O/c. This pseudomass ap- 
pears because the hamiltonian contains one large element 
at high energies, triggering a Lorentz-violating seesaw. 
Other models using combinations of mass and coefficients 
for Lorentz violation can be constructed to yield various 
exotic E^ dependences at particular energy scales. Note 
that the high-energy pseudomass and hence neutrino os- 
cillations depend on the declination Q of the propaga- 
tion. High-energy neutrinos propagating parallel to ce- 
lestial north or south experience the maximum pseudo- 
mass Amgo = a? /c, while others see a reduced value 
Atoq = Arngo cos^ 6. For propagation in the equatorial 
plane, all off-diagonal terms in /leff vanish and there is 
no oscillation. 

The features of atmospheric oscillations in the model 
are compatible with published observations. For defi- 
niteness, we take Amgo near the accepted range required 
in the usual analysis and Eq below the relevant ener- 
gies: Amgo = 10-3 eV2 and Eq = 0.1 GeV (c = 10-l^ 
a = 10-2° GeV). High-energy atmospheric neutrinos 
then exhibit the usual energy dependence, despite having 
zero mass differences. The zenith- angle dependence of 
the probability P^^^^^ averaged over the azimuthal angle 
also is comparable within existing experimental resolu- 
tion to a conventional maximal-mixing case with two gen- 
erations and a mass-squared difference Am^ = 2 x 10-^ 
eV^, as is shown in Fig. 1 for latitude x — 36°. However, 
the model predicts significant azimuthal dependence for 
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FIG. 1: -Pi^^^i/j, averaged over azimuthal angle for the bicycle 
model (solid) and for a conventional case with mass (dotted). 
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FIG. 3: (-Pi^5-.i^5 )adiabatic averaged over one year (solid) and 
at intervals between an equinox and a solstice (dashed). 
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FIG. 2: P^^_>^j, averaged over zenith angle for the bicycle 
model (solid) and for a conventional case with mass (dotted). 
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FIG. 4; (P^^_i,^)adiabatic for some modified models. 



atmospheric neutrinos, which is a signal for Lorentz vi- 
olation. For example, consider neutrinos propagating in 
the horizontal plane of the detector. Neutrinos originat- 
ing from the east or west have cos 9 = 0, Attiq = 0, and 
hence no oscillations. In contrast, those entering the de- 
tector from the north or south experience a pseudomass 
of Atoq = Ato§o cos^ X- Figure 2 shows the survival 
probability averaged over zenith angle as a function of az- 
imuthal angle. Although this model predicts no east- west 
asymmetry beyond the usual case, north-east or north- 
south asymmetries appear. Similar 'compass' asymme- 
tries are typical in all direction-dependent models. 

The basic features of solar-neutrino oscillations pre- 
dicted by the model are also compatible with observa- 
tion. Observed solar neutrinos propagate in the Earth's 
orbital plane, which lies at an angle r] ~ 23° relative to 
the equatorial plane. The value of cos^ Q therefore varies 
from zero at the two equinoxes to its maximum of sin^ 23° 
at the two solstices. Assuming adiabatic propagation in 
the Sun, the average Ve survival probability is 

(Piy^^iy, )adiabatic siu^ 9 siu^ Oq + COS^ 9 COS^ 6*0, (5) 

where is the mixing angle at the core, given by re- 
placing —cE with —cE + GpTie/V^ in Eq. (4). Figure 
3 shows the adiabatic probability as a function of en- 
ergy averaged over one year. The predicted neutrino 
flux is half the expected value at low energies and de- 
creases at higher energies, consistent with existing data. 
Also shown is the adiabatic probability at approximately 
weekly intervals between an equinox and a solstice. Over 
much of the year, it remains near the average. There is a 
strong reduction near each equinox, but the adiabatic ap- 
proximation fails there because oscillations cease, and so 
the true survival probability peaks sharply to unity. The 



combination of effects produces ripples in the binned flux 
near the equinoxes, which might be detected in detailed 
experimental analyses of existing or future data. 

Although detection of the semiannual variation would 
represent a definite positive signal for Lorentz violation, 
its absence cannot serve to eliminate this type of model. 
Simple modifications of the model exist that exhibit sim- 
ilar overall behavior for solar and atmospheric neutrinos 
but have only a small semiannual variation. As an illus- 
tration, consider the replacement of the coefficient (ciL)f^ 
with a coefficient (ol)^^ of half the size. This has the ef- 
fect of replacing the solid and dashed curves of Fig. 3 
with those shown in Fig. 4. The semiannual variations 
in this type of model lie below existing statistical sen- 
sitivities. Replacing also (a^)^^ with {aL)J^r is another 
option, which removes all orientation dependence in the 
model. Another example of a small modification is a 10% 
admixture of (ol) which raises the survival probability 
of 0.5 at low energies to about 0.6 without appreciably 
affecting other results. The ensuing survival probabil- 
ity in the adiabatic approximation is shown as the dot- 
ted line in Fig. 4. Other more complicated modifications 
that could be countenanced but that nonetheless retain 
the fiavor of the simple model include allowing depen- 
dence on directions other than Z, or even introducing 
arbitrary coefficients (aL)^^, {aL)'^^, {aL)tr^ and (cl)^|', 
which yields a model with 21 degrees of freedom. More 
general possibilities also exist [7] . We conclude that pos- 
itive signals for Lorentz violation could be obtained by 
detailed fitting of existing experimental data, but that it 
is challenging and perhaps even impossible at present to 
exclude the possibility that the observed neutrino oscil- 
lations are due to Lorentz and CPT violation rather than 
to mass differences. 
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The observations from long-basclinc experiments are 
also compatible with the oscillation lengths in the sim- 
ple two-coefficient model. For example, the oscillation 
length 27r/A3i controls survival and is short enough 
to affect KamLAND [26] . An analysis incorporating the 
relative locations of the detector and the individual re- 
actors would be of definite interest but lies outside our 
scope. Note, however, that the average survival prob- 
ability is {Pp^^pJ = 1 - 2 sin^ 6* cos^ 9 > 1/2. A complete 
analysis is therefore likely to yield a reduced flux some- 
what more than half the expected flux, in agreement with 
current data. 

The new class of long-baseline accelerator-based ex- 
periments [27] , planning searches for oscillations in j/^ at 
GeV energy scales and distances of hundreds of kilome- 
ters, will be sensitive to sidereal variations. The model 
predicts Uj- mixing with an experiment-dependent 



pscudomass Atoq = AmQo cos^ O because their beam- 
lines are in different directions and so involve a different 
propagation angle G. The energy dependence and tran- 
sitions will be similar to the usual mass case. 

Although the simple bicycle model reproduces most 
major features of observed neutrino behavior, it incorpo- 
rates only a tiny fraction of the many possibilities allowed 
in the SME. More complexity could be introduced in per- 
forming a detailed fit to all existing data. Nonetheless, 
the model serves to illustrate a few key phenomena in- 
troduced by Lorentz violation. It also shows that the 
presence of Planck-scale Lorentz and CPT violation in 
nature could well first be revealed by a definitive signal 
in neutrino oscillations. 
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